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Abstract 

Model matter agglomerations, with temperature as leading control parameter, have 
been considered, and some of their characteristics have been studied. The primary 
interest has been focused on the grain volume fluctuations, the magnitude of which 
readily differentiates between two commonly encountered types of matter agglom- 
eration/aggregation processes, observed roughly for high- and low-density matter 
organizations. The two distinguished types of matter arrangements have been de- 
scribed through the (entropic) potential driving the system. The impact of the 
potential type on the character of matter agglomeration has been studied, prefer- 
entially for (low-density) matter aggregation for which a logarithmic measure of its 
speed has been proposed. A common matter diffusion as well as mechanical relax- 
ation picture, emerging during the mature growing stage, has been drawn using 
a phenomenological line of argumentation. Applications, mostly towards obtaining 
soft agglomerates or so-called jammed materials, have been mentioned. 
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1 Introduction 



Agglomeration and aggregation of matter [1] are ubiquitous phenomena ex- 
pected to occur almost everywhere in nature. Their appearances range from 
metallic polycrystals to colloidal aggregates and clusters as well as lipid- 
protein viscoelastic matrices called often biomembranes. As observed in col- 
loids and protein aggregates immersed in solutions [2,3,4] and in granular 
media, such as (wet) sands [5], we can basically expect two types of emerging 
aggregations, namely, those termed diffusion-controlled, for which the attach- 
ment of particles and clusters occurs instantaneously at first touch, and the 
others, which are chemical-reaction controlled. For the latter, the 'first touch' 
sticking principle does not hold, and the microstructural evolutions have to 
follow certain attachment-detachment sequences of events. It finally results in 
expected matter incorporation in the form of a particle, or some cluster of 
particles, being attached to a colloidal fioccule or protein aggregate, for exam- 
ple. The resulting microstructures, however, are of different types [2,5]. Those, 
which are chemical-reaction controlled are rather more compact, whereas their 
diffusive counterparts look more irregular viz fractal, meaning that their arms 
or branches expand more visibly, so that they are quite undense as a whole 
[6,7]. A qualitatively similar observation can be offered while looking at matter 
agglomeration in some systems investigated in metallurgy [8]. Although one 
can find many qualitative similarities in the emerging structures, there are 
differences mostly due to different types of the matter fiuctuations, so that 
one could be able to distinguish among the emerging transformations [9,8]. 

We may easily observe that there are matter agglomerations that form typi- 
cally closely-packed (CP) microstructures. Polycrystals, ferroelectrics or ferro- 
magnets at sufficiently low temperatures, i. e. composed of electric or magnetic 
domains, respectively, fall into this category [10]. Binary alloys may also build 
structures in a domainwise way. In such systems, the total volume is a con- 
served quantity [11]. On the contrary, there exist processes, such as phase 
separations in solid solutions, as for example the sintering of powders [1], that 
lead to the formation of relatively loosely-packed (LP) microstructures such 
as immersion-containing assemblies, or certain composites [7]. These certainly 
do not occur under the constraint of constant total volume, see [12], and Ref. 
[25] therein. 

The potential governing the system, viz acting on the grains, plays an impor- 
tant role in a mesoscopic description. It has been shown for a type of model 
grain growth processes, that a logarithmic potential emerges, which is associ- 
ated with the pivotal role played by the curvature in such processes. The result 
obtained seems to be valid for many of the agglomerations whose emergence 
obeys the Laplace- Kelvin- Young law [13] . The results are also of importance 
for soap froths, foams and typical bubble-containing cellular microstructures. 
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droplets and microemulsions as well. All together they constitute the agglom- 
erates that we name CP microstnictures [14]. This seems to be in contrast 
to LP microstructures or aggregates. Based on the mesoscopic nonequilibrium 
thermodynamic description, in the present paper we are following an approach, 
that leads to a power form of the potential driving the system. This power 
form differs, however, from another genuine power form that we have derived 
previously and referred to as Lennard- Jones form throughout [15]. In fact, in 
an Euclidean space, in which the distance is defined by the grain radius, it is 
a linear form. Otherwise, it has algebraic form. 

The basic aim of the present study is to propose a model suitable to reveal 
essential differences between CP and LP types of matter agglomerates [14,12]. 
Since in [15,14], and in other studies [11,16], CP agglomerations have been 
studied quite exhaustively, in the present study we wish to place further em- 
phasis on LP agglomeration, describing it in detail. Moreover, we would like to 
offer a comprehensive picture of the matter agglomeration versus aggregation 
effects, understood roughly in terms of a sol-gel phase change, focusing the 
attention, however, not on the type of the phase transition but on the phys- 
ical factors controlling the phenomena in question. These are: temperature, 
space dimension as well as a type of interactions, but taken exclusively at the 
cluster level. To get an overview of the phenomena, we are also interested in 
viscoelastic properties of the slowly growing agglomerate in its final stage of 
evolution. 

A special motivation for the present study appears to be, that in the preceding 
studies [11,13,12,16] we derived the same growth rule for both CP and LP ag- 
glomerations. We found, however, certain differences between both processes: 
They do not behave in the same way as far as both, the number of grains as 
well as the total volume, are concerned [11]. Now, we wish to see which are the 
most essential differences between agglomerations grain volume fluctuations 
(T^{t). These show a different behavior because matter organization during 
agglomeration is different [17] under local temperature-damage conditions, as 
was also observed in fracture phenomena [18,19]. It seems to be almost exactly 
like in the thermal fuse model, designed for fracture phenomena, in which the 
same quantities, such as temperature and strain-stress factor play a leading 
role [19]. 

The plan of the paper is as follows. In Section 2, we propose a certain sum- 
mary of the principal characteristics of the two agglomerations/aggregations 
processes. We emphasize a simple derivation of a new (entropic) potential 
form and point to a similarity between two types of random walks. In Se- 
cion 3, we consider the grain volume fluctuations. In Secion 4, we discuss the 
overall context of the late-time growing phenomena of interest, coupled addi- 
tionally to mechanical stress relaxation, observed during late stage of growth. 
Conclusions are in Section 5. 
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2 CP and LP agglomeration processes and their characteristics 



In two recent studies [14,15] in which we have worked out a thermodynamic- 
kinetic description of model complex matter CP agglomerations, as well as in 
[12], where LP agglomerations (aggregations) have been considered, we have 
used the expression for the current in the space of cluster volumes [14] 



/, 



(1) 



where / = f{v, t) stands for the number density of grains of volume v, which 
means, that f{v,t)dv is a (relative) number of grains of a size taken from 
the narrow volume interval [v,v + dv] [13,11]; t is the time. represents the 
physical potential (see, [15] for an explanation of the term), assumed to be 
one of the relevant drivers for the agglomeration process at the mesoscopic 
level, and assuring its nonequilibrium character [14,15]. Moreover [14], 



K„) = ^„", . = 1-1, (2) 

where D is a diffusion reference constant; note that a principal role of D is 
to scale the time variable. kB represents the Boltzmann constant. Note that 
the mobility b{v) is related with the Onsager coefficient, L{v), that appears in 
derivation of the matter flux (1) under a set of assumptions, mostly on locality 
of the process, etc., for details see [14]. L{v), and hence b{v), could be measured 
by comparing the current and the thermodynamic force. The quantity Dv°' 
is to be inferred from the Green-Kubo formula, so that there would be quite 
a strong theoretical support for deriving b{v) both, experimentally as well as 
theoretically [14]. In [14,15] we have derived a (so-called) compaction potential 



4>{v) = 4>oln{v/vo), (3) 

where 0o, Vo - constants. While speaking about the role of the logarithmic 
potential, we wanted to see it as a quantity assuring emergence of rather com- 
pact structures [14,11]. Its role is equivalent to the action of grain curvatures 
[14] , resulting for a given value of the intergrain surface tension, in emergence 
of the driving force, such as the one involved in the Laplace-Kelvin- Young law 
[13]. 

In a previous study on the phase transformation kinetics for LP "difi:usive" 
agglomerates we have started with the matter flux of a purely diffusive nature 
[12,11], namely 
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(4) 



The diffusion function D{v) — Dv"" is proportional to the grain surface. 

Both CP and LP processes follow from the general form (1). Indeed, the case of 
LP is obtained when the second term in r.h.s. of (1) can be neglected. Formally, 
b{v) — > when T — > oo. Prom the physical point of view, it corresponds to 
sufficiently high temperatures T > T^h > 0, where Tfh can be treated as a 
threshold temperature above which the agglomeration takes place exclusively 
by yielding LP microstructures. However, the second term in r.h.s. of (1) 
depends both on v and T. It tends uniformly (i.e. independently of v) to zero 
at the high temperature limit if 



b{v)^^ const. (5) 
ov 

Then for a given system, temperature Tth does not depend on v and is consis- 
tently defined. It is a case when the potential 



4>{v) OC v^"" = v^/'^. (6) 

It should be noted that 4> is an entropic potential [15] which can be a cause of 
some disaggregation, or at least, cessation-to-growth (or, impingement) effects, 
taking place within the overall aggregation space. 

The above potential form (6), designed for LP agglomeration, seems convinc- 
ing as well as legitimate here: note that the 'force' 



rn-o OC — OC — OC -—r-r, (7) 

because v ~ , with R - the single grain radius. Thus, Fg^g acts as the in- 
verse of the area of the grain hypersurface which implies that the smaller the 
area is, the bigger the force acting on the grain can be, thus impeding the 
formation of new grains, which would contribute to an agglomerate's density 
increase. Note that this way we have established a certain matter densifica- 
tion procedure, expected to work in the systems of interest. Qualitatively, the 
similar dependence is for the CP matter aggregation: from (3) one gets for the 
'force' 



Here, Fg^g acts as the inverse of the hypervohime of the grain, which makes 
a difference between CP and LP agglomeration. Thinking in a somewhat dif- 
ferent way, some relaxation of the surface tension conditions [14] is expected 
to occur in the LP case, and some visible escape from the Laplace-Kelvin- 
Young strategy, c/. the review [11] for details, characteristic of polycrystals 
and cellular microstructures [13] can be foreseen. Moreover, the mean grain 
curvature may not be necessarily a dominating physical factor: either another 
type of curvature, the Gaussian curvature, may prevail [20], and the so-called 
model nano- agglomeration may occur [15]; or any relevant signatures of curva- 
ture contributions, as important thermodynamic-kinetic factors accelerating 
agglomeration, are not found at all [12], as occurs in some dispersive systems 
[3,21]. This is also what we mean by the relaxation of the surface tension 
conditions for loosely-packed grain systems. 

Referring further to (5) and using as in [15,11] the similarity relation, v ~ Ef-, 
one gets 



0(i?)oc0(i?o)^, (9) 



where 



(l>{Ro) = ^i?o, T > Tth, (10) 

and, in consequence, 4>{R) oc R; = D{1 — a), Rg as the initially observed 
grain radius is usually of the order of a few /im, or somewhat bigger [15,14], 
which imphes a typical length scale of colloids [3] . Notice, that while exploring 
the similarity relation ^ v ~ it!*^, one immediately sees that the state- dependent 
(grain-surface related) diffusion coefficient D{v) ~ D{R'^) oc R'^^^ is always 
obtained. Let us recall right here that for certain model biomembranes, that 
can be both reactive as well as constrained curvilinear systems, a modified 
kinetic coefficient, hke D{v, t) — Dv°-{t-\-tch)~^ , with tch a characteristic time, 
and h a microstructure-sensitive exponent [11], has also been proposed. The 
relation obtained can be anticipated as characteristic of systems with well- 
developed surfaces. 

^ As an extension ol v ^ R'^, one could think of clusters of fractal dimensions 
< df < d, aggregating in an Euclidean d-dimcnsional space: It is interesting to 
note that such cluster-cluster (colloidal) aggregation studies, with their structural 
as well as kinetic aspects involved, are mostly performed in d = 2, and for some 
of them, df ~ 1.6 was measured by several different methods [17]. Note that the 
potential (f>{R) from (9) will no longer be linear in R when a modification like 
V ~ R'^f ,0 < df < d would finally be used. On the contrary, it will take an algebraic 
(nonlinear) power form in R, with a characteristic exponent oi < df /d < 1. 
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Some supporting explanation of the surface-tension relaxation effect, from 
which the LP agglomeration may readily result, is under CP-agglomeration 
conditions. When the grains are packed densely together, the boundary of an 
individual grain can quite strictly be equivalent to a corresponding piece of 
the typically very thin overall intergrain separation space of the agglomerate. 
When the high-temperature activation conditions are switched on, the inter- 
grain separation space would tend to increase. This means that now the grain 
boundary is not as before, but very likely a second diffusive layer emerges, 
which effectively thickens its total width. Such a scenario is often associated 
with the double layer (DL) of Stern type, and bears a visible signature of 
phase separation effect [22]. 

Now, let us try out a phenomenology, supporting a plausible passage between 
the random walk in the space of cluster volumes [11,13], and the random 
walk, that is typically observed when studying the cluster-cluster phenomena, 
i.e. that of Einstein-Smoluchowski (ES) type [6]. In contrast to the agglom- 
eration/aggregation phenomena that we systematically propose to consider 
[12,14] in the space of grain/cluster sizes, the ES random walk is realized in 
the space of positions of the centers of mass of moving clusters [2]. For fur- 
ther remarks on a certain equivalence between both types random walks, and 
which are the (fractal) measures of it, see [16]. Why would one speak of such 
an equivalence ? Because both random walks in question can reasonably be 
anticipated as 'doing the same', at least in a thermodynamic sense: A (local) 
move of the grain or cluster boundary forward, associated with a step over 
an energetic barrier would be equivalent to an incorporation of the particle 
(macroion; molecular cluster), that, in turn, must overcome a thermal energy 
barrier to be a part of the boundary, or to be electrostatically caught by 
it. Otherwise, the boundary fluctuates, experiencing no advancement in the 
overall hypervolume space, so the particle also moves, performing its random 
walk. This is a typical scenario in a DL of Stern type and coincides well with 
cluster-cluster ES hke picture [22,23,16]. 

If we accept an ES-like picture of the cluster-cluster aggregation process, com- 
pleted by the involvement of DL, we can arrive at the following quantita- 
tive description while approaching the temperature threshold T^^. As can be 
learned from [24], where a complex viscoelastic behavior at the sol-gel transi- 
tion has been sketched as T ^ Tth from below, a sol phase, typically containing 
small clusters, the average linear sizes of which diverge at the mentioned tem- 
perature limit, merges rapidly into a bigger (gel) mega-cluster. Under such 
circumstances, according to [24] (see, Eq. (1) therein), the solvent viscosity, 
77s, involved in the Einstein-Stokes (E-S) relation, = ksT /Qt^tIsRh [25,22], 
where Rh stands for the hydrodynamic radius of the cluster, becomes cluster- 
radius dependent which explains the coagulation phenomena, often assisted 
by some changes in solution's viscosity [26]. Such a possibility can easily be 
inferred from our modeling as well. 
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Let us present it for spherical clusters. When applying forniTila (10) to the E-S 
formula, one is able, upon identifying Rh with Rg, to write 



Da = 11 

where, according to slow growth requirement postulated, \ R — Ro \ would al- 
ways tend to a small value for almost all t-s taken from the time zone of interest 
[13], thus enabling to make such an identification; it! = R{t). In consequence, 
an expression for the solvent viscosity can be proposed 



''O 



as well as the characteristic time t — l/rjs 



where Oq = AtkRo . 

By presenting the above picture we would like to point out a striking analogy 
between the sol-gel transition, and the LP-CP transition: Note that to accept 

the analogy, the temperature course has to be inverted here because, according 
to our previous explanations, the CP phase is usually a low temperature phase, 
and in consequence, LP emerges at higher temperatures. 



As was indicated in [24], the E-S diffusion coefficient depends on R^ but not on 
temperature, when the system is at the critical regime. The above described 
situation would invoke recent studies on the so-called kinetic arrest during 
gel formation, which can halt phase separation due to spinodal decomposition 
[27], thus enabling the gel formation. If the kinetic arrest is really the case, two 
physical quantities appear to be of relevance: the space dimension and the mag- 
nitude of the volume fraction wizthe aggregate density become important [28], 
apparently at the expense of domination of temperature, which is normally 
a prevailing physical factor [29,30]. A subtle interplay of attraction-repulsion 
fields [27], due to available electrostatics is a key phenomenon, so that our 
electrostatics-oriented picture of the percolative-in-nature process [28] is also 
supported this way. Such a distinction between reaction and non-reaction, viz 
diffusion limits is also very characteristic of colloidal (or, alike, such as those 
protein- involving [4]) aggregations, quite irrespective of the space dimension 
d [17,22,2]. 
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3 The grain volume fluctuations in the late-time growing stage of 
the agglomeration process 



There seems to be no doubt [8,7,3,2] that both agglomerations in question 
emerge in a fluctuating matter environment. Therefore, any reasonable quan- 
tiflcation of the so-called fluctuation impact on the speed of the matter evolu- 
tions under study is worth examining here. Thus, we propose below an eval- 
uation of the reduced variance 



as a direct measure of the grain volume fluctuations. 

The notation used in Eq. (14) refers to the statistical moments 



Kv^'it) > ^ j v''f{v,t)dv 71 = 0,1,2,... (15) 



of the agglomeration process, where the process is usually described by a local 
continuity equation 



dfiV:t) dJ{v,t) 


supplemented by the corresponding initial (of delta-Dirac type as a flrst suc- 
cessful attempt) and boundary (typically, of absorbing type) conditions (IBCs) 
[13,12,11,16]. The explicit solutions, /(v, t)-s, are presented elsewhere [12,13,11] 

Then zeroth moment, < v^{t) >, is related to the average number of grains 
in the system, and usually suffers from an algebraic drop with time [12,13]. 
The flrst moment, < v^{t) >, is related to its total volume which is a constant 
value for CP agglomerations [13,11] and an increasing function of time for LP 
agglomerations [12]. Prom the expressions of both moments, it follows that the 
average grain radius, Rav{t), behaves powerly with time, with growth exponent 
l/{d+ 1) that apparently carries some signature of random close-packing of 
matter by having included the super- dimension d + 1 [16,12,14,15]. These 
constitute the main characteristics of the model agglomeration/aggregation 
process in its late-stage {t » 1) limit. 

The question remains about asymptotic values of the moments < f "(t) > that 
must be known when applying formula (14). For CP agglomerations, they are 
found to obey a power law [13,11] 
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< V"{t) > ~ t^ri-l)/{2-a) = 0, 1, 2), t » 1, (17) 

whereas for LP agglomerations one provides another power law [12,11] 

< v"(i) > - i[("-i)+"]/(2-«) = 0, 1, 2), t » 1. (18) 

Notice, that for a = (d = 1) both power laws above converge to the 
same form, namely < v"(t) > ~ t("-i)/2_ When applying (14) and (17) it 
appears that for CP agglomerations, cr'^{t) can be fully identified with the 
inverse of < v^{t) > (the average number of grains), cf. [13,11] for details, 
what because of the constancy of < v^{t) >, leads to cr^(t) oc fsp(t), where 
(t) ~< v'^{t) > / < v^{t) >, and can be named the mean specific 
volume of the CP agglomerate, being equivalent to the inverse of its mean 
number density. The specific volume fluctuations, which are very important 
in nucleation phenomena [20,8], are 

a^{t)o^t''/^''+'\ (19) 
and if d — > oo, cr^(t) — Vgp oc t. 

When using (14) and (18), however, it turns out that for LP agglomerations 
(T^(t) is a quantity equivalent to the average grain radius Rav{t), see [12,11] 
too, and they behave in time as 

a\t)oct'/^''+'\ (20) 

and when d — > oo, (T^{t) ~ Ravit) — const. Note that a diffusional regime, 
characterized by the one-half exponent, is achieved exclusively for LP agglom- 
erations in ci = 1 because the only (linear) characteristic is Rav = Rav{t)'- Note 
that Vsp{t) is not a linear characteristic, since Vsp{t) oc [Rav{t)]^ usually holds, 
and here the d — 1-case must clearly be discarded as 'classically' diffusional, 
see Eq. (19) for comparison. 

Upon looking more carefully at the formula (19) and (20), we have found 
out an interesting 'harmonic', or efficiency, property for the agglomerations 
under study. Namely, we have observed, that upon defining a d-dependent 
logarithmic speed for the agglomeration in dimension d as 



we are able, up to a certain but sufficient accuracy level, to see that only for 
LP agglomerations one is capable of finding the harmonic rule, namely 
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^ +At, (22) 



V (2) 1/ (1) 1/ (3)' 
'^sp "^sp "^sp 



where simply, c/. Eq. (20) 



i^.p^'^ = ^, (i=l-3. (23) 

Thus, the case old — 2 looks as the most efficient from the aggregational point 
of view, since it contains mean (harmonic) properties of the d-dimensional 
process subject to its (speed) logarithmic fractal-like measure defined by (21). 
Many aggregation processes, sec [31,32], where colloidal monolayers are formed 
on the air-water interface, subjected to many external factors (UV light, AC 
nd DC currents, ions), and multitude of microstructures emerging in them 
(such as: fractal, bubble-like, mesomorphic, compact and voids-containing ag- 
gregates), have been experimentally observed and quantitatively analyzed, c/. 
[17], pointing to a Langmuir-Blodgctt-likc two-dimensional setup as the most 
frequently used device. For other interesting studies, here on the structure and 
phase transitions in Langmuir lipid monolayers, see a review [33]. 

To sum up in part, let us state that the fluctuations o"^(t) have been proposed 
as a reliable criterion of differentiating between CP and LP agglomerations 
in dimension d, with an emphasis placed on (i = 2, where harmonic-mean 
properties can be picked up in a natural way. It is known [34] that the di- 
mension d = 2, called a critical dimension, provides a natural framework 
for order-disorder phase changes in dipolar systems [22,23], where, in a low- 
temperature regime, one observes high-energy cost configurations, composed 
mostly of isolated dipoles, whereas at a high-temperature limit, a plasma 
(here, a low-energy cost configuration), consisting of some charged entities, 
emerges. Such a scenario can be compared to coagulation-fiocculation effects 
in colloids [3,31,32,17], so that the state of high T corresponds to flocculation 
phenomena - for them our approach with (9), taken at a mcsoscopic level [15], 
looks promising. Thus, a LP microstructure, of non-constant hypervolume, 
< v^{t) >7^ const [12], means a low-energy cost microstructure. Such a differ- 
entiation procedure can thus be anticipated as another mesoscopic procedure 
[29], helping in the selection of characteristic fluctuations in nonequilibrium 
instabilities and pattern formation [35], mostly those of hydrodynamic types. 



4 Matter diffusion and stress relaxation 



The above reasoning toward quantifying the fluctuations of the system can be 
strengthened with a phenomenological co-argumentation. The idea comes from 
a " coupled" diffusion-relaxation picture that appears in our complex system. 
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In any diffusion-migration growing process, the meclianical strain-stress fields 
play their role too. In our case, such a situation can be safely expected in the 
temperature domain T < T^^, though another type of relaxation of the stress 
field, say cr^, is expected to prevail when the CP agglomeration conditions 
are met. A different behavior may be observed when the CP agglomeration 
conditions are lost for the first time, that is, right at T = Tth, when the LP 
context appears. In both temperature regimes, the relaxation of crm{t) over the 
course of time, is very likely to go in a way, complementary to what the action 
of the matter diffusion-and-migration process, essentially described by the flux 
(1). This is expected to occur [14,12,11] presumably under homogeneous strain 
conditions, ~ const, for i >> 1. 

An additional motivation for drawing a coupled matter agglomeration and 
stress relaxation picture comes from the nature of the agglomeration process 
which can be considered to be the converse of its fracture counterpart. Con- 
sider a {am, Cm) plane on which the equilibrium curve is drawn, representing, 
for a given time instant, the state of a soft material under stress {am) that is 
elongated as measured by the strain (e^). After the quasi- linear Hooke regime, 
a nonlinear and plastic regime follows before failure. One can roughly consider 
the equilibrium curve to look like a Van der Waals (cubic equation of state) 
curve, meaning that the other extremity consists of the fragment region, which 
under pressure would agglomerate [19,18]. The region between the minimum 
and maximum on such a curve is known to be the siege of a first-order like 
transition region, characterized by nucleation and growth [18]. Moreover, the 
percolate nature of fracture process has been evidenced too, and its critical 
character, associated with the notion of phase transition of continuous type, 
with power law characteristics, has also been emphasized, mostly by numerical 
modeling of the process [19,10]. 

Prom [36] it can be learned, that in the absence of non-Arrhenius or fractal 
type kinetics, visibly modifying the diffusion coefficient D{v), one expects the 
Maxwell dashpot-and-spring model to represent the relaxation behavior. We 
wish to set up here a phenomenological picture, showing that both agglomer- 
ation and mechanical stress relaxation [37], where the stress relaxation takes 
place under slow growth conditions, proper of a mature growing stage, are cou- 
pled processes. To proceed, we will represent one of the two contiguous and 
matter-exchanging grains in the agglomerate, say grain 1, as an expanding one 
viz the spring, growing at the expense of its neighbor, to be termed grain 2 
(the dashpot), to which, according to the Maxwell model, the shrinking action 
should be assigned, , see [38] . 

For the system with narrow gaps, the Maxwell's model conditions are almost 
satisfied, so that the two-grain action can be extended over all pairs of contigu- 
ous grains until the expanding (growing) eventually survive. In a next step, the 
same kind of competition appears as in the well-argumented Laplace- Kelvin- 
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YoTing scenario suitable for cellular systems [13]. This picture holds in the CP 
context. 

In the LP context (a system with wide gaps), qualitatively we may have the 
same picture but with two differences which implies that grain expansion 
would not be likely so vigorous that the corresponding gap is wider, therefore 
untight, so that the fluid leakage might a priori be more pronounced. Thus, 
the fluid response against the piston wall is weaker, and the stress relaxation 
of Maxwell type no longer applies. 

The stress relaxation can now be described by introducing a relaxation expo- 
nent X ill the Maxwell-like, or quasi-fractional Maxwell [38], model presented 
below. This exponent is, in general, d-dependent. Coupled matter diffusion 
and stress relaxation processes have been analyzed in [39] for an anomalous 
random walk in a position space. Here, we wish to offer a similar picture, but 
for our random walk performed in the space of cluster volumes. As is known, 
the Maxwell stress relaxation picture leads to an exponential decay of the 
stress: am{t) ~ exp^—t/r), where r is a reference time of the concentrated 
grain-containing viscous solution, Eq. (13). This behavior holds for T < Tfh- 
As mentioned previously, for T > Tfh we propose 



^ + ^ = 0, (24) 
at T 

where the above is usually true when the internal strain field, is practically 
constant, which is the situation expected in rheological model behavior of 
colloids and polymers, see [40], and refs. therein. When solving (24), one arrives 
at 



am{t) {t/ry'/^''-'\ t»l, (25) 

where x = 2d -|- 3 = 2{d -|- 1) -|- 1 ^ The obtained formula, Eq. (25), resembles 
to a large extent an empirical formula cr^(t) oc t~'^, see Eq. (9.30) in [41], 
characteristic of mechanically relaxing amorphous polymers; a dimensionless 
exponent m is a measure of the relative importance of elastic and viscous 
contributions, in our model manifested exclusively by a close relation with d, 
a rough messure of losses in the system cf. Eq. (25). Usually m is proportional 
to the loss tangent, a well-known dissipation factor in such processes [41], so 



^ It is interesting to note here that x = x{d) can be, upon identifying d as con- 
secutive numbers of (pairwise) emerging bonds in a gelling system, a certain odd- 
number based generator of the Bethc lattice elements, starting from the 3-bond 
(initial) generator for d = 0, and continuing with d, what perfectly recovers the 
number of elements of model Bethe lattice (an ultrametric space) , a tool very useful 
for mean-field description of gels, and other bonds-containing systems [22,5,2]. 
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that the energy dissipation is due to the character of the agglomeration (and 
relaxation) space, characterized by x involved in Eq. (25). 



The overall relaxation exponent in (25) reads 



—1 = wrry P*" 

which is exactly one half of the growth exponent I'sp^'^'' given in (20), sec (23) 
too. Such a result can be exactly obtained under the assumption of vaUdity 
of the Hall-Petch (HP) relation [42]. For T Tth [36], one recovers a volume 
fluctuation-driven growth-and-relaxation synchronous mode 



CTm ~ i?a.-'/' ~ (27) 

which clearly confirms the solution presented above, Eq. (25), cf. relation 
(20). Notice that HP relation typically holds for low-temperature systems [42], 
though its range of validity still expands because of its applications to new 
studied systems [42,36]. Here, in fact, we conjecture a validity of the HP rela- 
tion outside the low-temperature range, for example in (dense) colloidal gels, 
or other (reactive) polycrystalline and/or jammed materials (zirconia, that 
means, ceramics produced by sol-gel method), in which sol-gel characteristics 
are very much related to polycrystalline behavior [43,44]. 

Although such a conjectured supposition needs to be proved, it is found to 
play a significant role in many, mostly plastic as well as superplastic (also, an 
inverse HP relation, characteristic of superplastics) metallurgical systems [8]. 
It also holds for non-metallurgical systems, such as (bio) polymeric for example, 
cf. [36,11], and refs. therein. Proceeding this way, the HP relation couples, at 
a local thermodynamic equilibrium point near Tth, diffusion and migration 
of matter, as well as a complex material relaxation under slow growth in a 
fluctuating medium, especially in its mature stage, t >> 1. 

Another relaxation exponent can be defined 



where 



/^./^ = 2(dTl) ' ^=1'2'3> (29) 
holds, and a very simple relation, like 
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"^I^sp^ ^ — ^sp^ \ d — 1) 2, 3, 



(30) 



comes out, in fact, as a consequence of the HP condition applied for our 
complex electrorheological context. 

Similarly, the harmonic rules favoring the dimension d = 2, in which most of 
interesting experiments arc performed [33], also apply when looking at /isp'-'^-', 
so that the relaxation in dimension two appears to be an expected case and 
well worth experimenting. Coated materials, recently produced by sol-gel tech- 
nique deserve recognition, cf. [44], and refs. therein. In consequence, it is be- 
lieved, that such a complex diffusion-relaxation behavior may readily help in 
explaining many clustering phenomena taking place under carefully designed 
experimental conditions, see [31,32,17] as selected cases. In such experiments, 
an account of the stress micro-fields, arising in a direct relation with elec- 
trostatic fields acting on all the charged objects cannot be underestimated, 
either at short or even at long separation distances, see [45]. In [31,32] such 
a micromechanical factor, being a cause of differentiation amongst the ob- 
tained latex-sphere- containing microstructures, and the possible influence on 
fractal or non-fractal properties of the resulting mega-clusters, often leading 
to different coverage values, has not yet been taken into account. 



5 Conclusions 

The present study is devoted to a comprehensive picture of the model matter 
agglomeration/aggregation phenomena at a mesoscopic (cluster) level. The 
integrated picture has been achieved by offering, mostly in a phenomenological 
way, some characteristics that are: temperature, space dimension as well as 
interaction-type dependence, and completed by some mechanistic point of 
view [37], due to coupled-to-agglomeration relaxation of matter during late- 
time growing stage. This can undoubtedly be understood as a promising new 
contribution of our study. 

The grain volume fluctuations, generally given by Eq. (14), and in particular, 
by Eqs. (19) and (20), for CP and LP agglomerations, respectively, appear 
to be a reliable criterion for understanding the fluctuations in two different 
types of model systems, conforming, however, by construction [12,14,11], to 
the same average grain radius involving growth rule, namely the one given by 
Eq. (20). It proves then that the growth rule is not a sufficient criterion to 
differentiate between two seemingly different types of processes, even in their 
mature stages. 

A logarithmic potential (3) (p{R) oc ln{R/Ro) [15], can also be derived based 
on an analysis of the electric component of the interfacial bending moment 
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presented in [46] for micelles (cylindrical molecular condensers). A power- 
like Lennard- Jones potential, obtained in [15], also follows from that analysis, 
and has been derived for micelles of another type, called spherical condensers 
[46,40]. In the present work, we have picked up a potential 0, Eq. (9), in its 
simplest form, i. e. for non- fractal objects, being linear in R (and, recognized in 
elementary electrostatics as well, but with a clearly specified viz. temperature- 
as well as dimension-dependent prefactor, (j){Ro)), Eq. (10), certainly being of 
entropic character [29]. 

It is interesting to note that the harmonic-mean rule, Eq. (22), applied for 
the logarithmic speed defined by (21), has been detected: It seems that any 
matter compacting CP "procedure" would somehow destroy certain efficiency 
conditions, and a natural appearance of spontaneously formed mesostructures. 
This may also be the case when not only d-dependent dimensional constraints 
apply but also when the system is subjected to external physical factors that 
may change the number of the degrees of freedom, causing certain stress-strain 
fields induced losses in the coverage value, cf. [31,32], and the pictures therein. 

An overall experimentally-motivated [1,2,4,8,37] picture that emerges from 
our model study shows that the diffusion-migration process may thoroughly 
be accompanied by mechanical relaxation, and that both mutually cooperat- 
ing sub-processes may go synchronously. Note that any approval for such a 
synchronous way implies that must also ffuctuate, because Rav is identified 
with cr^(t) for a LP process, so that the quasi-fractional Maxwell model pro- 
posed (24) is likely to work under the requirement that this quantity fiuctuates 
"accordingly", cf. Eq. (27), i.e. the overall synchronisation mode. In conse- 
quence, a simple relation, interconnecting both sub-processes described, may 
emerge in such a simple form, see (30). The harmonic- mean rule becomes appli- 
cable in both coupled cases. Both processes are recognized as inseparable (re- 
call sol-gel coexisting phases [44,24]), and go in the same direction; notice, the 
"-" sign in (28), which is, however, not the case of (21). Invoking here a com- 
parison presented in the previous section between agglomeration/compaction 
and fracture/compaction processes (being mutually converse), we see also cer- 
tain well-established parallels between our model and the scalar anti-phase 
mechanical model of rupture with elastic interaction, cf. [19,18]. 
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